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1. Introduction 

This article represents the second part of our series on cosmological models within 
alternative gravity theories. In the first part M, we presented a cosmological model 
in Weyl-Cartan spacetime which is no longer tied to a Riemannian spacetime and is 
based on an alternative gravity theory, so-called metric-affine gravity (MAG), which 
was essentially developed by Hehl et al and has been extensively reviewed in ||. In 
|| we derived the field equations of our new model and provided exact solutions for a 
rather broad class of parameters. We were able to show that it is possible to construct 
models in which the non-Riemannian quantities die out with time. 

Within this work we investigate the observational consequences of the new model, 
and try to rule out if it is compatible with current observations. As becomes clear 
from the title we are interested in the derivation of the so-called magnitude-redshift 
relation, which establishes a relation between the luminosity and distance of a given 
object and the parameters within the underlying cosmological model. This relation was 
successfully used to extract parameters within the cosmological standard model, also 
known as Friedman-Lemaitre- Robertson- Walker (FLRW) model |l|, ||, which is based 
on General Relativity (GR). 



Recent observations O, O, 13, 14, 15] of type la supernovae, which are believed to 



be some kind of standard candles, led to the unexpected discovery that we seem to live 
in an universe which is currently undergoing an accelerated phase of expansion |T3L [T5 



The availability of a data set incorporating over 90 type la supernovae has turned the 
magnitude-redshift relation into a cosmological test which must be passed by any new 
model. It interesting to test if predictions, like the one that the expansion of the universe 
is currently accelerating, still hold in case of our Weyl-Cartan model. Therefore we put 
special emphasis on the model dependence of certain assumptions which are made during 
the derivation of the magnitude-redshift relation. It is noteworthy that the analysis 
of the SNe data led to a renewed interest in the cosmological constant A within the 
cosmological standard model. Since the cosmological constant seems to be a concept 
which re-emerges periodically within cosmology it is interesting to figure out if it is also 
inevitable for the description of the data within an alternative model. Note that there 
are also several efforts, cf. p7|, |29], Q e.g., to cope with the new observational situation, 



which do not focus on a change of the underlying cosmological model. 

Although parts of this work are based on the results obtained in Q, it should also 
be of use for the reader who wants to learn more about the derivation and the use of 
the magnitude-redshift relation within the standard FLRW model. At this point we 
like to mention that the model in Weyl-Cartan spacetime still has toy model character, 
since many important issues are not yet worked out. Hence this article represents a first 
attempt to rule out whether our new model is viable when confronted with the currently 
available observational data from the type la supernovae. 

The plan of the paper is as follows. In section we review some standard concepts 
of the FLRW model and introduce the notions and assumptions which are necessary to 



A cosmological model in Weyl-Cartan spacetime II 



3 



derive the magnitude-redshift relation within this model. In section [3] we make use of 
these results in order to derive this relation within the new Weyl-Cartan model. We 
perform fits within the FLRW as well as within the Weyl-Cartan model to the combined 
data set of Wang |24| in section f|. Finally, we compare our results with the results of 



several other groups who performed similar analysis within the FLRW as well as in 
non-standard model and draw our conclusion in section [5]. In Appendix A , we provide 



an overview over the units used throughout the preceding sections. 

2. Magnitude-redshift relation within the cosmological standard model 



In this section we present a short derivation of the magnitude-redshift relation within 
the standard FLRW model. The results obtained here will be of use when we perform 
our analysis within the new Weyl-Cartan model in the next section. In particular we will 
stress the model dependence of the magnitude-redshift relation in the following section, 
especially its entanglement with the field equations. Since the upcoming derivation 
of the field equations can be found in every standard textbook on cosmology we keep 
ourselves rather short. For further reading the reader is referred to QH [I]. 



2.1. Field equations 

The assumption of homogeneity and isotropy leads to the so-called Robertson- Walker 
metric as starting point. Using spherical coordinates (t, r, 6, (ft) and the coframe 

= dt, $ l = dr, $ 2 = S(t) r d9, $ = S(t) r sin 9 d<f>, (1) 

V 1 — kr 2 

the line element is given by 

ds 2 = -0 a ® 0° + tf 1 <g> $ l + d 2 <g> d 2 + <g> s . (2) 

The function S(t) is the cosmic scale factor, and k can be chosen to be +1, —1, or 
for spaces of constant positive, negative, or zero spatial curvature, respectively. The 
only thing missing in order to set up the field equations is an appropriate matter model. 
Following the cosmological standard model, we assume that matter is smoothly smeared 
out over the whole spacetime. Thus, we choose the energy-momentum 3-form of an ideal 
fluid with pressure p, energy- density /i, and four- velocity u a , i.e. 

sr=w, (3) 

where r] a := *d a and 

Z a P = ( fJ , + p)u a uP +pg af3 . (4) 
The Einstein field equations with cosmological constant read: 

r! af3l A R*> + 2\r] a = 2 K £« uid . (5) 

Here R a p denotes the Riemannian curvature 2-form, r] a/31 := * (0 a A ^ A 7 ), A is the 
cosmological constant, and k is the gravitational coupling constant. Insertion of the 
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Roberston- Walker metric yields the following set of field equations (the dot denotes 
differentiation with respect to t) 

= g/^ (6) 

k 

— - A = - Kp. (7) 

Hence the field equations (H) turned into a set of ordinary differential equations for the 
scale factor S(t). The set of additional unknown parameters A, k, fi, and p depends on 
the model we decide to consider. Note that \i and p are related by an equation of state 
(eos) p = 

In addition to the field equations we have one Noether identity, which takes 

the following form in a Riemannian spacetime 

DH a = 0. (8) 

Let us assume that the equation of state takes the form p(t) = w n (t), with w = const. 
Using ([1|)-(FJ), and (|3|), equation (||) turns into 

fiS = -3S (/I + p) P ^ ii = X!S- 3(1+U,) . (9) 

Thus, we have found a relation between the energy density and the scale factor, which 
depends on the constant w in the equation of state and an integration constant K\. 

2.2. Critical density 

In case of a flat universe, k = 0, the Hubble function H := S/S and the density n are 
related via a unique function [i c = n(H), which is often called the critical density. The 
critical density is obtained via the first Friedman equation (|Bp, in case of a vanishing 
cosmological constant we havef| 

u2 k k 8ttG k=0 3H 2 

H + lP=r = —» * ^^^G- (10) 
This quantity is called critical density because of its character to distinguish between a 
open, flat, or closed universe, i.e. 

^ = { R2 + '^j^G ^ Pk=-i < < Pk=i- (11) 

As we can infer from ([10]) the critical density is determined by measuring the current 
value of the Hubble function. Let us introduce the dimensionless density parameter Q w , 
which represents the relationship between the actual and the critical density 
fi 8tiG k 

nw = ]T c = MP lX= (12) 



| Note that we make use of natural units (cf. Appendix A ), i.e. h = c = 1. Thus, the gravitational 
coupling constant becomes k — 8irG. 
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Note that we use w as an index since we did not specify the underlying equation of state. 
In case of a non-vanishing cosmological constant A the first Friedman equation might 
be written in terms of a total density parameter which encompasses both contributions, 
i.e. fi total = tt w + tt x := 3§?/i + 3^, yielding 

R2 + ^-^ = ^ ^ tal -l = ^. (13) 

Note that the value of total density parameter distinguishes between the three possible 
geometries of the 3-dimensional subspace, i.e. 

{< 1 ( < open 

= 1 ^kl =0 flat . (14) 
> 1 \ > closed 

2.3. Redshift 

Maybe one of the most striking properties of the cosmological standard model is the 
redshift due to the global expansion. In the following we will focus on its derivation by 
using our ansatz for the metric from equations (|I])-(0). 

The first assumption we make is that the propagation of light coming from a distant 
object, a galaxy e.g., can be treated as a classical wave phenomenon. With r = r\ being 
the radial coordinate of the object, r = r coordinate of the observer, t — t\ time 
at which wave is emitted, and t — to at which the light is detected at r = r , the 
Robert son- Walker line element yields 

, 9 „ #=<A=const f dt f 1 dv . . 

i8 =° L wrL tttw- (15) 

Thus, we gained a relation between the coordinate distance and the time, setting r = 
we obtain 

p dr f t0 _dt_ 

Jo VT^kF 2 ~ nn) ~ k S(ty 

Shifting the time of the emission t\ —>■ t\ -\-5t\ will result in a shift of the detection time, 
i.e. to — > to + <^o- Since the lhs of (|i"6| ) does not depend on the emission or detection 
time we can infer 

to dt f t 0+ 5t ^ rt +St ^ rh+St! ^ 

^ I 777TT = / TvTV (!7) 



fiTi) = / ~7u7\- (16) 



Sit) J tl+Stl S(t) J t0 S(t) J tl S(t) 



St ,i « t 0A -> S(t) = const = (18) 

Interpreting 5t ,i as t ne times between two wave crests at emission and at detection, i.e. 
relating them to the wavelength <5t ,i ^ Ao,i, one obtains 
Ai S(h) 



A s{t y 

Using the common astronomical definition of the redshift, i.e. 



(19) 



Ao 1 J > redshifted ,„„x 
~ ~* 1 < blueshifted ' ^ } 
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we find the following relation between the scale factor and the redshift in a Robertson- 
Walker spacetime 

1+ *=m (21) 

Thus, if the universe is expanding, then distant sources should be redshifted. Note that 
the derivation presented above makes use of three crucial assumptions: (i) the paths of 
photons are completely determined by the metric, (ii) the dispersion relation \v = c is 
valid at all times, and (iii) the observed wavelengths are small compared to the size of 
the universe. 

Let us now derive an alternative form of the Friedman equation. For that purpose 
we make use of the expression for the redshift (pTj) , and the relation between the scale 
factor and the density @. With Q k : = —^57^2 equation fll3|) turns into 

tt w + tt x + tt k = 1. (22) 

For normal matter (w = 0) we can rewrite the Friedman equation as follows 
2 k A (|) 2 K k A 



} } H 2 = Hi 



i> 0-2 



r/j m0 (1 + zf - fc (1 + zf + 



H 2 = Hl[(l + z) 2 (l + zn m0 ) -z{2 + z) n xo ]. (23) 

Here and throughout the rest of the paper we denote present values of certain quantities, 
i.e. at the time to, by an index 0. The Hubble parameter in terms of the redshift reads 

H = Jt l0g ( So J = Jt bg ( IT! J = —zdt (24) 



fj ^ _ u-i ,1 r,i - ,2 

dz 



Hq 1 (1 + z)- 1 [(1 + zf (1 +zOno) -z(2 + z) n xo ] 5 . (25) 



Before we proceed with the derivation of the luminosity distance in the FLRW model, 
we note that there is a remarkable connection between the density parameters and the 
so-called deceleration parameter q. The deceleration parameter, which is commonly 
introduced when one expands the scale factor around a certain time (cf. [|IJ), is defined 
as follows 

S S S d 1 . . 

»== - = -JFs = Jt H ~ L < 26 » 

Thus, in case of a FLRW model which contains only usual matter and a cosmological 
constant, cf. equation (p3|) , the deceleration factor at present time is given by the simple 
expression 

g = % - Q xo . (27) 
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2-4- Luminosity distance 

In order to assign a distance to objects one introduces the so-called luminosity distance, 
using the fact that the light from objects far away from us appears fainter than the light 
from nearby ones. The astronomical definition reads 

energy per time produced by source \ 2 
energy per time per area detected by observer J 

( luminosity \ 1 _ ( L \ 2 
47r x flux J \ AttF ) 



^luminosity 



{ 



Thus, by measuring F and with knowledge of L (via a standard candle, a supernova 
e.g.) we are able to determine the distance <i luminosity . Of course this distance definition 
implies that we know how much light is emitted by the source at least during a specific 
time interval. We will not investigate this question here any further since the search for 
an appropriate model of the source belongs to the realm of astrophysics (a discussion of 
type la supernova models can be found in [p6fl ). The question we have to ask ourselves 



is: How is the luminosity distance related to parameters within the FLRW model? Since 
energy is conserved the following equation is supposed to hold 

L 5tx 5Xi = F 5t 5X A , (29) 

here A denotes the area of the 2-sphere at the detection time t = to, and Sto,i, an d ^-Vi 
the different length and time scales at emission and detection due to global expansion. 
The distance d FhRW between an object at r = r%, which emits light at t = ti, and an 
observer at r = r = 0, who detects the light at t = to, is given by 

d^ RW = S{t )r 1 Ao = 4nd 2 FhRW = 47rS 2 (t )rl (30) 

Thus, with the help of equations flT9|)-(]2"TD the observed flux can be expressed in the 
following form 

f = - 5 Jl 5Xl (19) ±© - n + 7]- 2 ) 

4irS 2 (to)r 2 5t 5\o ~ A7rS 2 (t )r 2{ 1 ' K ' 

Comparison with (128|) yields 

rfluminosity = S (t ) T\ (l + Z) . (32) 

Additionally one wants to replace r\ by the scale factor S. We make use of equation 
(|16|) and obtain 



/(n) = / jf= = I *° -^y = { n for k = o 




arcsin (ri) k — +1 

r\ for k = 

arcsinh(r!) k = — 1 



k = +1 

ri=< T^T for k = 

k = -l 
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sm 



d 



luminosity 



S(t ) {l + z)x 



*o dt \ 
tx S(t) I 



(33) 



k = +1 
for ^ = • 

Thus, we are able to express r\ via the scale factor. In order to derive r± explicitly we 
need a solution of the Friedman equations, which of course depends on the cosmological 
model, i.e. the choice of the parameters in (§|)-(0), we decide to consider. Note that we 
did not make use of the field equations of the underlying gravity theory up to this point. 
This fact will be crucial when we derive an expression for the luminosity distance within 
a cosmological model which is not based on general relativity in the next section. As we 
will show in the following subsection there is an elegant way to rewrite the luminosity 
distance in terms of the density parameters, which makes use of the expression for the 



Hubble parameter as derived in equation (|23"D . Note that this is the point when the field 
equations come into play. 

Special case In case of a Friedman model which contains only normal matter and a 
contribution from the cosmological constant one can express the luminosity distance as 
a function redshift and the model parameters. As we will show below, this expression 
is used when one wants to perform fits to observational data. Again we make use of the 
Robert son- Walker line element, i.e. 

dr dt (Eh So 



VI - kr 2 S 



VI - kr 2 



dr 



dr 



o \/l — kr 2 



;i + z) dt, 

dz 



dz 



z (2 + z) tt 



AO 



i + zf(i + Z Q m0 )-z(2 + z) n xo 



d 



luminosity 



s (i + z) e 



[HqSo)- 1 [ F[z\d~z 
Jo 

Where we made use of the following definitions Q[x] : = 



(34) 



sin (x) k — +1 

x for k = and 

sinh(x) k = — 1 

F[z] := [(1 + z) 2 (1 + z Q m o) — z (2 + z) ^ao] 2 • If we make use of the definition of f4 
we can rewrite equation (R3) as follows: 

(1 + *) 



dr 



luminosity 



-^ovl^fcol 

(1 + z 



|fifco| / F[z]dz 
Jo 

Vll-^O-^Aol 



e 



F \z) dz 



.(35) 



Ho\/\l — f2 m o — ^ao| 

Thus, within a Friedman model with normal matter and cosmological constant the 
luminosity distance turns out to be function of the corresponding density parameters, 
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Table 1. Assumptions made up to this point. 



Ansatz /Assumption 



Equation 



General Relativity as underlying gravity theory (5) 
Metric is of Robertson- Walker type (i.e. homogeneity and isotropy) (I)-(^) 

Photons follow null curves (i.e. are determined by the metric) (15) 

Observed wavelengths small compared to size of the universe (IS) 
Dispersion relation valid at all times, and peculiar movement 

of the source neglect able (1£) 

Sources of known constant absolute magnitude ( 36 ) 
Photons travel unimpeded between source and observer, i.e. no 

gravitational potentials or dust between source and observer (@)~@ 



the Hubble constant, and the redshift, i.e. <i luminosity = d luminosity (z, H , f2 m0 , Q\o). Note 
that this is a remarkable result since G?i uminosi t y depends only on the present day values of 
the parameters within the model and the redshift. 

2.5. Magnitude-redshift relation 

Due to historical reasons astrophysicists often use the so-called magnitude as unit for 
the luminosity of a stellar object. The relation between the distance-redshift relation 
and the so-called magnitude-redshift relation (cf. ^ |TT|, 0) is given by 

m(z, H , Q , Q xo , w, M) := M + 5 log (^^) + 25 

\ length / 

= M + 5 log (H d lnminosity ) - 5 log f-^-] + 25. (36) 



^ length t 

Where M represents the absolute magnitude of the observed star. By introducing a new 
constant M. := M — 5 log if o + 25 we are able to express the distance-redshift relation 
in equation (|3~5"D in a compact way as magnitude-redshift relation 

m(z, H , O m0 , ttxo, M) 



M -i 5 leg I ; + ^ =9 

v/|i - n m0 -n A0 | 



Vll-OnO-flAol / F [z] dz 
Jo 



■ (37) 



This relation is commonly used to extract cosmological parameters, like the density 
parameters associated with normal matter and the cosmological constant, by performing 
fits to data sets which were produced by the observation of standard candles, i.e. objects 
of known absolute magnitude. Equation fl3T|) will be of use in section |] where we will 
perform fits to a real data set. Note that table [I] contains a collection of all assumptions 
made during the derivation of the magnitude-redshift relation within the cosmological 
standard model. 
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3. Magnitude-redshift relation within the Weyl-Cartan model 

After our review of the derivation of the magnitude-redshift relation within the 
cosmological standard model we will now switch to the new Weyl-Cartan model which 
was presented in [[|. We will not discuss the details of this model at this point, therefore 
the reader should consult sections 2-4 of |§ in order to get an idea which assumptions 
were made during the derivation of the field equations of this model. 



3.1. Field equations 

Let us collect the field equations derived in || (cf. eqs. (56)-(62)). Note that we make 
use of the form in which the constant S, which was introduced in eq. (54) of [Q, is 
already set to zero. 




(a 4 + a 6 )K — 




S 

+ (a 4 + a 6 ) k | | — 




(38) 




(39) 

A, (40) 

(41) 




Additionally, we have the following relation between the pressure and the energy-density 
fi = 3p r -8c(^] . (42) 



As one can see from equations (j38|)-(|4~2|) there is an additional function ((t) entering the 
field equations, besides of the scale factor S(t) from the Robertson- Walker metric. This 
function stems from an ansatz for the Weyl 1-form Q which governs the non-Riemannian 
features of the model, cf. section 4 of ||. Here k denotes the usual parameter within 
the Robertson- Walker metric, x> a 4? a 6, c, and n are coupling constants, and A is the 
so-called induced cosmological constant. 
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3.2. Magnitude-redshift relation 

As we have shown in sections [2.4| and |2.5| , the magnitude-redshift relation depends on 
the mechanism of light propagation of the underlying theory of gravity and the field 
equations. Since we now work in a space which carries nonmetricity, besides of the 
torsion which does not affect light propagation, the question arises if the trajectories of 
photons deviate from the ones in the Riemannian theory. The answer to this question 
is no, since in a Weyl-Cartan spacetime the null geodesies are determined by the usual 
geodesic equation as in the Riemannian case. Hence we are allowed to use relation 
([2~4|), which is a consequence of the Robert son- Walker metric, in order to derive the 
magnitude-redshift relation. Thus, our next aim is to express the Hubble function in 
terms of some density functions and model parameters. From equation (f4~0D we infer 



H 2 



S - - 

's's 2 ' 



In order to eliminate the second order term we make use of equation 
the solution for S/S from this equation into (|43| ) leads to 



H 2 



(A 2 S 2 - 2Ak) k (a 4 + a 6 ) - x& + nS 2 p r - 4kc( 5 



(2kA (a 4 + a 6 ) + x) S 2 



(43) 

I). Reinserting 
(44) 



H 2 
11 o 



(2kA (a 4 + a 6 ) + x) 

2 



Xk 



tt2 o2 



Sq 
3 




(45) 



(0) 



HI 



(2kA (a 4 + a 6 ) + x) 



-X^ho (1 +z; + T7 2 



k (04 + a 6 ) [HqCI 
AkcC, 2 



2 

AO 



2Afi fc0 (1 + zf] 



~7H TF202 \\ + z ) 



tt2 02 
n o 



In the last equation we introduced the density parameters 



H 2 S 2 , 



and Q? 



(46) 



A 

H 2- 



Subsequently, we have to choose an equation of state and an appropriate ansatz for (. 
We choose the eos to be of the form p = wfi with w = const, and make use of the 
solution for jj, obtained in || equation (65), i.e. n = —jzr§^ (§) ■ Hence we can infer 



H 2 



Hi 



(2kA (a 4 + a 6 ) + x) 



k (a 4 + a 6 ) [Hln 2 m - 2AQ k0 (1 + zf] 



fcO 



1 



w 



4kc( 2 

H 2 S$ \l-3w 



(47) 



For ( we make use of the solution mentioned in 
const, which finally yields 



equation (70), i.e. ( = i/S with l 



H 2 



Hi 



(2kH 2 Q A0 (a 4 + a 6 ) + x) 



- X n k0 (1 + zf - 4fi C0 (l + z 



(a 4 + a 6 ) [HlQ 

1 — w 
l-3w 



2 

AO 



2-^o^AO^fcO 



(48) 
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where we introduced the new density parameter fir '■ 



C ■ H 2 S 4 ' 



Special case In case of a vanishing induced cosmological constant, i.e. A = 0, equation 
(filf) turns into 



X 



4(1 + z ) 4 n 



Co 



w 



1 



X^lko (1 + z) 



(49) 



1 -3w / 

In order to calculate the magnitude-redshift relation we remember equation (p4|) and 
insert this expression into the first line of equation (|34]) , yielding the following luminosity 
distance 



^luminosity^, H Q , flkO, tyo, X; W ) 



:i+z) 



in 



AO | 



X Jo 



G[z]dz 



with 6 as defined after eq. fl3§, and G[z] := [4 (1 + if fi ?0 (f^) - X ^w (1 + zf] 
Hence the magnitude-redshift relation is now given by (cf. eq. (|5SD) 

m(z, # 0) ^co ; ^ko, w, M, x) = M + 5 log | H d luininosity (z, H Q t0 , tt k0 , x, tu)} 



(50) 



.M + 5 log 



e 



\Q 



fcO | 



(51) 



vT^oi |_V x Jo 

Note that in contrast to the FLRW model there is no simple relation as in ( ^2) between 
the density parameters. Therefore we cannot eliminate the density parameter f4 in the 
equation for the magnitude. In case of a flat model equation (|51]) reduces to 

m(z,HQ,VL^w } M } x) 



M + 5 log 





7* 




</ 



4(l + z) 4 U 



CO 





i 

2 






j 


VI -3w/_ 







(52) 



In the next section we will make use of equation (^Tj) in order to determine whether it 
is possible to describe the available type la supernova data within our model. Finally, 
we mention that the Hubble parameter found in equation fl49D enables us to derive the 
deceleration parameter, which was defined in equation (|26|) . Its present day value is 
given by 

4fi (0 (w - 1) 



4. Numerical results 



(53) 



In this section we will present the numerical results obtained by fitting the magnitude- 
redshift relations in (|37| ) and ( j52]) to a real data set. We start with a collection of the 
different available data sets of type la supernovae, which were also used by other teams 
to determine the cosmological parameters. 
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Table 2. SNIa data sets. 



Symbol Number of SN Reference 



Comments 



I 

II 
III 

IV 

V 

VI 



18 
42 
10 
10 
1 

27 



p. 571, |1C,|11| 

p. 570, M 

p. 1021,111 

p. 1020, Jl|] 

@ 

p. 1035, H 



Calan/Tololo survey 
Supernova Cosmology Project 
High-z Supernova Search Team 
Same as III but MLCS method 
Farthest SNIa observed to date 
Low-redshift MLCS/templatc 
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4-1. Data sets 

In table ^| we collected the number of supernovae and the references which actually 
contain the data. Note that the data sets of the different groups is not directly 



comparable. Perlmutter et al [|TTJ provide the effective magnitude 



in the B band, 



while Riess et al [15] use the so-called distance modulus [M As shown by Wang in [24 



it is possible to find a relation between this two data sets by comparing the data of 18 
SNe la published by both groups. The definition of the magnitude as given in equation 
(j36p is compatible with the definition used by Perlmutter et al , it is related to the 
definition of Riess et al by 

m = M + /i = M + 51og d luminosity + 25 = M + 5 log H d luminoslty . (54) 

As shown in \UM we have to choose M = — 19.33±0.25 in order to transform the different 



data sets into each other. Note that this value corresponds to the MLCS method of Riess 
et al . In the following we will make use of the data of Wang which contains 92 data 
points and can be viewed as a compilation of the sets I, II, IV, and VI from table |2] in 
which some outliers were removed. 



4.2. Fitting method 

Since we want results which are comparable to the analysis of the combined data set by 
Wang in P3j| , we are going to minimize]]}] 



X 



92 



ll eory (zi [parameters) — y ( i ™ easurcd 



a z ■ + a A ■ 

(ii 1 mz % 



(55) 



in order to obtain the best-fit parameters within the standard and the Weyl-Cartan 
model f| H 0- Here /i* heory denotes the distance modulus at a certain redshift z, t as 
defined in fl5~4]) . The error of the measured ^™ easured is given by a 2 i . The dispersion 
in the distance modulus a mz due to the dispersion in the galaxy redshift, a z , can be 
calculated iteratively by 



ZralO 



1 



dd 



luminosity 



_ ^luminosity 



dz 



(56) 



Not to be confused with the energy-density within the field equations. 

Not to be confused with the coupling constant within the field equations in eqs. 
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Table 3. Grids used for minimization. 



Parameters 


[Interval, Stepsize] 




{f2 m o, Q\o} 

{n k0 ,n CQ ,x, w} 


{[-2. ..4, 0.01], [-2. 
{[-2. ..4, 0.01], [-1. 


.4,0.01]} 

.1,0.01], 1,[-1. ..0,0.1]} 



according to Wang (cf. equation (13) of ||24j| ). We perform a brute force calculation on 
grids as denoted in table |3] in order to find the minimum of (|55|). 



4-3. Best- fit parameters 

In table f| and |5] we collected the best-fit parameters obtained by the method described in 
the previous section. Note that we did not impose any constraints on our parameters, 
like spatial flatness, e.g., when performing our search. In figure [I] and |2] plotted the 
corresponding distance modulus versus redshift relation together with the data set of 
Wang which consists out of 92 type la supernovae. Note that the ^-distributions 
displayed in figure |] and |2] correspond to the plane containing the best-fit parameters 
found on our grid. In both figures the 95.4% confidence level corresponds to the outer 
boundary. In table |] and [7] we collected some of the results of other groups which 
performed a similar analysis within the FLRW as well as in non-standard cosmological 
models. Note that these collections are by no means exhaustive, we therefore apologize 
for not having mentioned all the works which are devoted to this subject. 



5. Conclusions 



Fitting results As we have shown in the previous section it is possible to describe 
the observational data within both models. It is noteworthy that without any 
additional constraints the best-fit parameters, i.e. set F13 in table |], within the 
FLRW model favour a closed universe, whereas the best-fit solution within the 
Weyl-Cartan model, i.e. set C2 in table corresponds to an open universe. In 



case of the FLRW model our results comply with the ones of Wang found in [24 



If one wants to impose the condition of spatial flatness, the best-fit parameters 
are given by (F28, 65, 0.29, 0.71, 135.26, 1.51, -0.56) within the FLRW model, and 
(C9, 69, 0, -0.28, 1, -0.5, 292.821, 3.32, 1) within the Weyl-Cartan model0]. This is an 
interesting result since as soon as we assume that the universe is flat we are not able 
to find parameters within our new model which fit the data. This is in contrast 
to the FLRW model where the assumption of spatial flatness worsens the fit only 
slightly. Hence, at least within the parameter intervals we considered, our model does 
not support a flat universe. Note that the fit F28 is in compliance with the current 
cosmological concordance model which encompasses about 30% of matter and a dark 
energy contribution of about 70%. 

^[ Note that here we make use of the same enumeration as in table [| and table ||. 
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Table 4. Best- fit parameters (FLRW model). 



uj I11UU1 


^0 


n 


1 'AO 


X 




X 


2 

V 




Fl 


50 


0.45 


-2.00 


401. 


.04 


4 


.50 


2.22 


F2 


54 


-0.60 


-2.00 


256. 


.81 


2 


.88 


1.70 


F3 


55 


-0.81 


-2.00 


230. 


.86 


2 


.59 


1.59 


F4 


56 


-1.00 


-2.00 


208. 


73 


2 


.34 


1.50 


F5 


57 


-1.18 


-2.00 


190. 


.31 


2 


.13 


1.41 


F6 


58 


-1.34 


-2.00 


175. 


.50 


1 


.97 


1.33 


F7 


59 


-2.00 


-2.00 


188. 


.04 


2 


.11 


1.00 


F8 


60 


-1.15 


-1.49 


155. 


.83 


1 


.75 


0.91 


F9 


61 


-0.64 


-0.79 


148. 


.80 


1 


.67 


0.47 


F10 


62 


-1.14 


-1.21 


145. 


.58 


1 


.63 


0.64 


Fll 


63 


1.02 


1.29 


142. 


,53 


1 


.60 


-0.78 


F12 


64 


3.10 


3.35 


233. 


.03 


2 


.61 


-1.80 


F13 


65 


0.63 


1.10 


134.54 


1.51 


-0.78 


F14 


66 


0.80 


1.40 


134. 


,91 


1 


.51 


-1.00 


F15 


67 


0.92 


1.63 


137. 


,02 


1 


.53 


-1.17 


F16 


68 


0.99 


1.80 


140. 


,98 


1 


.58 


-1.30 


F17 


69 


1.04 


1.94 


146. 


,84 


1 


.64 


-1.42 


F18 


70 


1.07 


2.05 


154. 


66 


1 


.73 


-1.51 


F19 


72 


1.06 


2.17 


176. 


,32 


1 


.98 


-1.64 


F20 


73 


1.04 


2.20 


190. 


.23 


2 


.13 


-1.68 


F21 


74 


1.02 


2.23 


206. 


,18 


2 


.31 


-1.72 


F22 


75 


0.98 


2.23 


224. 


.21 


2 


.51 


-1.74 


F23 


76 


0.95 


2.24 


244. 


,30 


2 


.74 


-1.76 


F24 


77 


0.91 


2.23 


266. 


,47 


2 


.99 


-1.77 


F25 


78 


0.87 


2.22 


290. 


,71 


3 


.26 


-1.78 


F26 


79 


3.05 


3.68 


733 


,20 


8 


.23 


-2.15 


F27 


80 


0.13 


1.45 


374. 


.68 


4 


.20 


-1.38 



[H ] = kms _1 Mpc _1 . 



Table 5. Best- fit parameters (Weyl-Cartan model). 



Symbol 


H 


f^/c0 




X 


w 


x 2 


xl 




CI 


65 


-1.07 


0.05 


1.00 


-0.80 


138.034 


1.56 


-0.109 


C2 


66 


-1.03 


0.05 


1.00 


-0.90 


138.028 


1.56 


-0.110 


C3 


67 


-1.01 


0.03 


1.00 


-0.10 


138.034 


1.56 


-0.111 


C4 


69 


-0.95 


0.04 


1.00 


-0.60 


138.056 


1.56 


-0.106 


C5 


70 


-0.92 


0.04 


1.00 


-0.70 


138.034 


1.56 


-0.105 


C6 


71 


-0.90 


0.03 


1.00 


-0.20 


138.031 


1.56 


-0.111 


C7 


74 


-0.83 


0.04 


1.00 


-0.80 


138.039 


1.56 


-0.113 


C8 


80 


-0.71 


0.03 


1.00 


-0.50 


138.035 


1.56 


-0.112 



[Ho] = kms _1 Mpc _1 . 
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Figure 1. ^-distribution in the (f2 m o, r2Ao)-plane and magnitude-redshift relation in 
case of the standard model. On the top lhs the x 2 -distribution in the H — 65 plane is 
shown, which also contains the best-fit parameter set F13. On the top rhs we plotted 
the distance modulus versus the redshift in case of the best-fit parameter set. The 
green dots correspond to the experimental data for 92 type la SNe as contained in the 
data set of Wang. In the two other figures we plotted the distance modulus versus the 
redshift for the other parameter sets in table 0. 



-0,9 -0.95 -1 -i.05 -i.i -1.15 -1.2 -1.25 



%0 




H = 66 




Figure 2. ^-distribution in the (QkO> £\o) _ pl ane and magnitude-redshift relation in 
case of the Weyl-Cartan model. On the lhs the ^-distribution in the Ho — 66 plane 
is shown, which contains the best-fit parameter set C2. On the rhs we plotted the 
distance modulus versus the redshift in case of all parameter sets in table ||. Since all 
of this sets fit the data very well the eight different curves appear as a single line at 
the selected resolution. 
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Table 6. Best- fit parameters other groups (FLRW model). 



Symbol Ref. Best-fit parameters 



Comment 



SI 

S2 

PI 

Rl 

R2 

Gl 

G2 

G3 

VI 

V2 

V15 

Wl 



[13| 
[13| 

[15| 
[15| 
[23| 
[23| 
[23| 
[19| 
[19| 
[21| 
[24| 



{n mQ = -o.2, n XQ - o} 

{O m0 = 0.4, n xo = 0.6} 
{O m0 = 0.28, n XQ = 0.72} 
{O m0 = 0.24, O A0 = 0.72} 
{O m0 = 0.2, n xo = 0.8} 

{o m0 = -o.i, n xa = o} 

{O m0 = 0.4, n xo = 0.6} 

{O m0 = 0.3, n xo = 0.7} 

{O m0 = 0.28, fl XQ = 0.72, M = 23.94} 

{O m0 = 0.79, n XQ = 1.41, M = 23.91} 

{O m0 = 0.33, n xo - 0} 

{#0 = 65,0^ = 0.7,0*0 = 1-2} 



MLCS 
Template 

MLCS 
Template 



1997ff included 
Combined data set 



[H ] = kms^Mpc -1 . 



Table 7. Best- fit parameters other groups (non-standard models). 



Symbol Ref. Best-fit parameters 



Comment 



Bl 

V3 

V4 

V5 

V6 

V7 

V8 

V9 

V10 

Vll 

V12 

V13 

V14 

V16 

Tl 



[22| 
[19| 
[19| 
[19| 
[19| 
[19| 
[19| 
[20| 
[20| 
[21| 
[21| 
[21| 
[21| 
[21| 
[25| 



{O m o 
{O m o 

{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{O m o 
{Oao 
{*! = 

Fo = 64,0^ = 0.6, Xtf = 0.3} 



0.3, O Rigid0 = 0.7} 
0.49, n xo = 0.51, Al = 23.97} 
1.86, n xo = 1.52, M = 23.95} 
0.4, Oao = 0.6, X = 23.96} 
0.98, n xo = 1.53, M = 23.91} 
0.4, Oao = 0.6, M = 23.96} 
1.62, O\o = 1-59, M = 23.93} 
= 0.46, M = 24.03} 
= 1.34, M = 24.03} 
= 1.41,^ = -1} 
= 1.22, w = -1} 
= 0.48} 
0.4} 
« = 5} 

0.87, O^ = 0.3, 



Conformal model 



= 0.54, Oao 
= 1-76, Oao 
= 0.79, O^o 
= 0.65, fi^o 
= 0.52, Oao 
= 0.6, Oxo = 
= -0.358, z ma 
0.08, H^/Hl 



X 
A 
A 
A 
A 
A 



5- 2 
5- 2 
H 2 
H 2 
M 
A* 

Variable A 
Variable A 
Quintessence model 
Quintess. model +1997ff 
A - H 2 +1997ff 
A ~ S~ 2 +1997ff 
QSSC model 
Model with local void 



[H ] = kms -1 Mpcf 



Deceleration factor The best-fit sets F13 and C2 in both models predict a universe 
which is presently in an accelerating phase of expansion. In case of the Weyl-Cartan 
model the current value of the deceleration parameter is roughly seven times smaller than 
in the FLRW case, hence in our model the expansion of the universe seems to accelerate 
less rapid than predicted by the standard scenario. Figure f| provides an overview over 
the sign of the deceleration factor in the FLRW and in the Weyl-Cartan model. In case 
of the latter we plotted the distribution in the parameter plane which contains the best- 
fit C2. Another interesting property associated with deceleration factor is the fact that 
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it seems to provide an independent cosmological test. As displayed in the plot on top 
lhs in figure |3| the contour lines for a constant deceleration factor within the (fi m0 , ^ao) _ 
parameter plane of the FLRW model intersect the ones for constant flko at a non-zero 
angle. Hence if we assume that we can measure Q^o via analysis of the fluctuations 
within the CMB then we are able to pin down the pair (f2 m o,f2Ao)- The situation 
within the new model is similar to the one FLRW case, i.e. the curves of constant go an d 
flko intersect each other at a non-zero angle. Figure ^| contains plots for five different 
choices of the equation of state parameter w. As we can see from plot at the bottom 
lhs only the choice w = 1 corresponds to a degenerated situation, since in this case the 
deceleration factor vanishes according to (|53|) . 



Other groups As we have shown in sections |2| and ^| the magnitude-redshift relation 
depends on several assumptions. Especially its strong dependence on the underlying field 
equations renders it to be a useful tool to discriminate between different cosmological 
models. In table || and [F] we collected some of the results from several other groups who 
used this relation within the FLRW as well as in non-standard models. It becomes clear 
from table |^ that in addition to the model dependence of the magnitude-redshift relation, 
the estimates for the cosmological parameters strongly depend on the data set which is 
used to perform a fit. It is interesting that some of the early best-fit parameter sets like 
SI and Gl, e.g., correspond to unphysical models. The situation within non-standard 
models is similar. As one can infer from table ^ the parameters estimated within non- 
standard scenarios depend strongly on the underlying model. Although it seems to be 
possible to describe the supernova data equally well within several alternative scenarios, 
the main benefit of the cosmological standard model consists of its level of detail, which 



results in the availability of a number of different cosmological tests [Tfj, |17|, |18| . Note 
that without the constraints from other cosmological tests the used data set does not 
favour a flat universe. The best-fit F13 is even slightly better than the one which is 
obtained if we impose the constraint of spatial flatness right from the beginning, i.e. 
the fit F28 mentioned at the beginning of this section. It is noteworthy that nearly all 
groups come to the same result that a non- vanishing cosmological constant seems to be 
inevitable for the description of the supernova data within FLRW model, as becomes 
clear from table ^| and our results in table |j. 



Summary & Outlook In figure [5] we plotted the distance modulus versus the redshift 
up to z = 2. The upper curve corresponds to the best-fit within the Weyl-Cartan model. 
It becomes clear that the supernovae at high redshifts will appear dimmer within this 
model than in the FLRW case. The data point at z — 1.7 corresponds to the farthest 
known supernova 1997ff as reported by Riess et al [14]. Although not significant this 
supernova seems to favour the best-fit FLRW model F13. The uncertainties connected 
with 1997ff are large, therefore we did not include it in our fitting procedure. Note 
that also the best-fit flat model F28 seems not to fit the data from 1997ff. A possible 
magnification of 1997ff by gravitational lensing was discussed by Mortsell et al in |28[] . 
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Figure 3. Contour lines of the deceleration factor qo and the density parameter 
fijto m the (f2^ j^mo) an d (^fccb^co) plane, respectively. The figure on the top left 
corresponds to the standard model which contains only usual matter w = and a 
contribution from the cosmological constant. The other figures belong to the Weyl- 
Cartan model in case of different choices of the equation of state parameter w and a 
vanishing induced cosmological constant. 

It is too early to make a reliable prediction at this point, one has to wait until more 
data at high redshifts becomes available. As we can see from the rhs in figure^ a survey 
at high redshifts (like SNAP e.g.) should enable us to discriminate between the two 
best-fit models C2 and F13. 

It is also possible to derive the age of the universe within our model. In case of 
the parameter set C2 we obtain an age of 18.4 Gyr. Thus, our model combined with 
the supernova data yields an older universe than the standard model - i.e. 13.8 Gyr in 
case of F13, and 14.6 Gyr for the flat scenario F28. Hence, there is no conflict with the 
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Figure 5. On the lhs we plotted the distance modulus versus the redshift for the best- 
fit models F13, F28, and C2 up to z = 2. As one can see from the rhs the difference 
between both models grows with the redshift. 



estimates from the observations of globular clusters and nucleocosmochronology which 
range from 11 to 15 Gyr |31, |32 . 



We reviewed the derivation of the magnitude-redshift relation within the 
cosmological standard model with special emphasis on the various assumptions which 
are necessary to obtain this relation. We have shown that it is possible to find a 
similar relation between the luminosity and distance within the new Weyl-Cartan model 
proposed in ||. We performed fits to the combined data set of Wang [24] within both 
models. Note that this is the first time that a fit of MAG based model to a real data 
set has been performed. Our efforts can be viewed as a first step to relief the Weyl- 
Cartan model of its current toy model character. We are aware that at the moment 
our model is far from being called realistic. Nevertheless we were able to show that in 
principle it is possible to obtain a magnitude-redshift relation within our MAG based 
model. Despite this success it is too early to use the best-fit parameters obtained in 
this work in order to pin down the free parameters within the model. One needs at 
least another independent cosmological test in order to impose meaningful constraints. 
Therefore, our aim in the next article of this series is to check whether our model is 
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Table Al. Natural units. 



[energy] [mass] 


[time] [length] 




length -1 length" 


-1 length length 




Table A2. Units of quantities. 




Quantities 


1 


Length 


Gauge potentials 
Field strengths 
Matter currents 


[Qap], [Ra/3] 


P.]- 1 


Coordinates 


[0]Mdr] 


M 

L J 


Functions 


[n w ],[n k },[n m },[si x ], 
[C], [fitotai], [fid. [°a] 




Miscellany 


[m], [M], [M], [u] 


4 5 [^luminosity] 5 


Constants 


M,M,[fc],[a/] 


M^[A]-i[A]-5,[ c ]-i, 

[S]-4,[G]2,[^ 1 ]-4 + 3(l + ™ ) 



compatible with recent observations of the cosmic microwave background and if there 
exists a concordance region with the parameters obtained in this work. 
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Appendix A. Units 

In this work we made use of natural units, i.e. H = c = 1 (cf table |A1|) . Additionally, we 
have to be careful with the coupling constants and the coordinates within the coframe. 
In order to keep things as clear as possible, we provide a list of the quantities emerging 
throughout all sections in table |A2| . Note that [d] = 1 and [*] = length n ~ 2p , where n = 
dimension of the spacetime, p = degree of the differential form on which * acts. 
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